Self-similar Magnetoresistance of Fibonacci Ultra-Thin Magnetic Films 
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We study numerically the magnetic properties (magnetization and magnetoresistance) of ultra- 
thin magnetic films (Fe/Cr) grown following the Fibonacci sequence. We use a phenomenological 
model which includes Zeeman, cubic anisotropy, bilinear and biquadratic exchange energies. Our 
physical parameters are based on experimental data recently reported, which contain biquadratic 
exchange coupling with magnitude comparable to the bilinear exchange coupling. When biquadratic 
exchange coupling is sufficiently large a striking self-similar pattern emerges. 

PACS numbers: 75.70.Cn; 75.70.Pa; 75.70.-i; 71.55.Jv 



The discovery of quasicrystals in 1984 [y aroused a 
great interest, both theoretically and experimentally, in 
quasiperiodic systems. One of the most important rea- 
son for that is because they can be defined as an inter- 
mediate state between an ordered crystal (their defini- 
tion and construction follow purely deterministic rules) 
and a disordered solid (many of their physical properties 
exhibit an erratic- like appearance) On the theoreti- 
cal side, a wide variety of particles, namely electrons 
phonons [Q, plasmon-polaritons spin waves etc, 
have been studied. A quite complex fractal energy spec- 
trum, which can be considered as their basic signature, is 
a common feature of these systems. On the experimental 
side, the procedure to grow quasiperiodic superlattices 
became standard after Merlin a/ [Q , who reported the 
realization of the first quasiperiodic superlattice follow- 
ing the Fibonacci sequence by means of molecular beam 
epitaxy (MBE). 

Parallel to these developments in the field of quasicrys- 
tals, the properties of magnetic exchange interactions 
between ferromagnetic films separated by non-magnetic 
spacers have been also widely investigated jsj. The dis- 
covery of physical properties like antiferromagnetic ex- 
change coupling dl , giant magnetoresistance (GMR) , 
oscillatory behavior of the exchange coupling |]ll[, and 
biquadratic exchange coupling (EEC) ||l2|, made these 
films excellent options for technological applications and 
attractive objects of research. For example, GMR in 
magnetic films has been widely considered for applica- 
tions in information storage technology . 

It is known that GMR also occurs in non-periodic gran- 
ular systems, like Cu-Co alloys consisting of ultrafine Co- 
rich precipitate particles in Cu-rich matrix p^ . Due to 
the fact that precipitate particles of these heterogeneous 
alloys have an average diameter and an average spacing 
similar to magnetic films, the origin of GMR in gran- 
ular systems is also similar to the one found in mag- 
netic films p5[. Therefore, quasiperiodic systems which 



present magnetoresistive properties can be a first step 
for a better understanding of magnetoresistance in gran- 
ular systems. On the other hand, from a technological 
point of view (as we will show later in this letter) the 
BEG associated with quasiperiodicity permit us to con- 
trol magnetic field regions, where magnetoresistance re- 
mains almost constant before saturation. 

The aim of this work is to investigate the influence 
of quasiperiodicity on the magnetic properties of ultra- 
thin magnetic films. In particular, we are interested 
in Fe/Cr(100) structures, which follow a Fibonacci se- 
quence, whose experimental magnetic parameters were 
recently reported by Rezende et al p^ . 

A Fibonacci structure can be grown experimentally by 
juxtaposing two building blocks A and B following a Fi- 
bonacci sequence. In our specific case we choose Fe as 
the building block A and Gr as the building block B. 
A Fibonacci sequence S'at is generated by appending the 
N — 2 sequence to the TV — 1 one, i.e., Sn = Sn-iSn-2 
{N > 2). This construction algorithm requires initial 
conditions which are chosen to he Sq — B and 5*1 = A. 
The Fibonacci generations are: 5*0 = [B], Si = [A], 
S2 = [AB], S3 = [ABA], etc. Therefore, the well known 
trilayer Fe/Gr/Fe is the magnetic counterpart of the third 
Fibonacci generation (A/B/A). We remark that only 
odd Fibonacci generations have a magnetic counterpart, 
because they start and finish with an A (Fe) building 
block. Fig. 1 shows schematically the third and fifth 
Fibonacci generations and their magnetic counterpart, 
where t {d) is the thickness of a single Fe layer (single Gr 
layer). It is important to note a double Fe layer whose 
thickness is 2t in the fifth generation corresponding to 
a double letter A. It is easy to show that the quasiperi- 
odic magnetic films, for any Fibonacci generation, will be 
composed by single Cr layers, single Fe layers and double 
Fe layers. 

We consider the ferromagnetic films with magnetiza- 
tion in the plane xy and take the z axis as the growth 
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direction (see Fig. 1). The very strong demagnetization 
field, generates by tipping the magnetization out of plane, 
will suppress any tendency for the magnetization to tilt 
out of plane. The global behavior of the system is well 
described by a simple theory in terms of the magnetic 
energy per unit area i.e., 

Et = Ez + Ebl + Ebq + Ea, (1) 

where Ez is the Zeeman energy, Eu is the bilinear en- 
ergy, Ebq is the biquadratic energy and Ea is the cubic 
anisotropy energy. It is usual to write the total magnetic 
energy in terms of experimental parameters (or effective 
fields) of each interaction, 

F " 1 
irf = V(t.A){-i?ocos(0, - Oh) + -HaSui\2ei)} 
tMs ^ 8 

n-l 

+ {-Hu cos{0, - 9,+i) + Hbqcos^ie, - (2) 

i=l 

where t is the thickness of a single Fe layer and we assume 
Mi = Ms. Hbi is the conventional bilinear exchange 
coupling field which favors antiferromagnetic alignment 
(ferromagnetic alignment) if negative (positive). We are 
concerned here to the case Hu < because magnetore- 
sistive effects occur only for this case. Hbq is the BEC 
field, which is responsible for a 90° alignment between 
two adjacent magnetizations and is experimentally found 
to be positive . Ha is the cubic anisotropy field which 
renders the (100) direction an easy direction. Hq is the 
external in-plane magnetic field and 9h is its angular 
orientation. From now on we consider 9h — 0, which 
means that the magnetic field is applied along the easy 
axis. The thickness and the angular orientation of the 
i-th Fe layer are given by U and 0i, respectively. 

The equilibrium positions of the magnetizations {9i} 
are numerically calculated by minimizing the magnetic 
energy given by (2). It should be remarked that it has 
proved difficult for us to generate accurately configu- 
rations for larger structures, mainly when the BEC is 
strong [ prf . However, we got results in sufficiently large 
generations to infer important informations about the ef- 
fect of the quasiperiodicity. 

Theoretically, the spin-dependent scattering is ac- 
cepted as responsible for the GMR effect |l^. It has 
been shown that GMR varies linearly with cos(A6') when 
electrons form a free-electrons gas, i.e., there is no barri- 
ers between adjacent films |l^. Here, A9 is the angular 
difference between adjacent magnetizations. In metalic 
systems like Fe/Cr this angular dependence is valid and 
once the set {9i} is found, we obtain normalized values 
for magnetoresistance, i.e, 

ri-l 

R{Ho)/R{0) = 5] [1 - '^^^(^^ " ^^+l)]mn - 1), (3) 

i=l 



where R{0) is the resistance at zero field. 

Now we present numerical calculations for the magne- 
tization and the magnetoresistance curves for Fibonacci 
ultra-thin magnetic films. The physical motivation for 
that is because the Fibonacci quasiperiodic structure can 
exhibit magnetic properties not found in the periodic 
case and the long range correlations induced by the 
construction of this system are expected to be reflected 
someway in the magnetoresistance curves. We have con- 
sidered physical parameters based on realistic values of 
the magnetic coupling fields, whose whose experimental 
data were recently reported |l^. We assume the cubic 
anisotropy field Ha = 0.5 kOe, corresponding to Fe(lOO) 
with t > 30 A growth by sputtering . We choose the bi- 
linear and the biquadratic fields, Hbi and Hbq, such that 
their values lie in three regions of interest: i) close to the 
region of first antiferromagnetic- ferromagnetic transition 
where Hbi is moderate ||l^; ii) near to the maximum of 
first antiferromagnetic peak, where Hbi reach its maxi- 
mum value 0; and iii) in the second antiferromagnetic 
peak, where Hbi is small and equal to Hbq . 

In Fig. 2 we show the curves of the normalized magne- 
tization and magnetoresistance versus the magnetic field, 
for the third Fibonacci generation (corresponding to the 
Fe/Cr/Fe trilayer). We assumed Hbi = —0.15 kOe and 
Hbq — 0.05 kOe {\Hbi\ > Hbq). These parameters corre- 
spond to a realistic sample with Cr thickness equal to 15 
A. From there one can identify two first order phase tran- 
sitions at Hi ~ 100 Oe and H2 ~ 220 Oe. Also there are 
three magnetic phases presented: i) an antiferromagnetic 
phase {Ho < 100 Oe); ii) a 90° phase (100 Oe < Hq < 220 
Oe); ui) and a saturated phase {Hq > 220 Oe). We re- 
mark that our numerical calculations indicate that a first 
order phase transition occurs when Ha > 2{\Hbi\ + 2Hbq). 
Since the transition magnetic fields are the same for both 
the magnetization and the magnetoresistance, from now 
on we concentrate our discussion on the magnetoresis- 
tance curves, because it is easier to investigate their self- 
similar pattern. 

Fig. 3 shows the normalized magnetoresistance curves 
for the fifth(a) and seventh(b) Fibonacci generations with 
the same experimental parameters considered in Fig. 2. 
In Fig. 3(a) we can identify four first order phase transi- 
tions, where each one is due to a 90° jump of magnetiza- 
tion. This behavior is always displayed when the BEC is 
present in the magnetic energy. Previous works on phase 
diagrams have looked carefully the origin and features of 
the so-called 90° phase [ p7|JlS| ]. For the seventh genera- 
tion there are eight first order phase transitions and nine 
magnetic phases are present from the antiferromagnetic 
phase {Hq < 38 Oe) to the saturated one {Hq > 440 
Oe). Note a clear self-similar pattern of magnetoresis- 
tance curves by comparing Fig. 2 and Fig. 3, i.e., the 
pattern of the trilayer Fe/Cr/Fe is always present in the 
next generations. On the contrary, when Hbi = — 1.0 kOe 
and Hbq = 0.1 kOe {\Hbi\ » Hbq), which correspond 
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to a sample with Cr thickness equal to 10 A, the self- 
similarity is not observed, as it is shown in Fig. 4. For this 
set of parameters, the majority of phase transitions are of 
second order and we have found numerically that this oc- 
curs when < 2{\Hu\ + '2Ht,q). However, when the ratio 
between Hbq and Hbi is increased {\Hbi \ = H^q = 35 Oe), 
we observe again a striking self-similar pattern (see Fig. 
5), where each new transition occurs for a value of mag- 
netic field which is about a half of the previous one. For 
this set of parameters the magnetoresistance is approxi- 
mately 1/2 its value at zero magnetic field, because the 
magnetizations of the adjacent Fe films are nearly per- 
pendicular to each other due to the strong biquadratic 
field. For the third generation. Fig. 5(a), there is only 
one transition at Hi ^ 70 Oe and two magnetic phases: 
a 90° phase at Ho < 70 Oe and a saturated phase at 
Tfo > 70 Oe 111. In the fifth generation. Fig. 5(b), there 
are two transitions at Hi ^ 70 Oe and H2 ^ 140 Oe, 
respectively. For the seventh generation, as one can see 
from Fig. 5(c), there are three transitions at Hi ~ 35 Oe, 
H2 - 70 Oe and H3 - 140 Oe. 

From the numerical results above discussed, we can 
infer that the magnetoresistance exhibits a self-similar 
behavior when: (a) Hbq is comparable to Hti, and (b) 
there is a first order phase transition (see Fig. 3 and Fig. 
5). A possible explanation for that is because the BEG 
reinforces the quasiperiodic order, which is responsible 
by the self-similarity in quasiperiodic systems. This is 
an unexpected effect of this unusual exchange coupling 
and, as far as we know, this is the first system which 
presents magnetoresistance with self-similar properties. 
Besides, from a technological point of view, magnetore- 
sistance with almost constant regions (Fig. 3 and Fig. 5) 
opens new perspectives in information storage technology 
by the possibility of a recording system with more than 
two states. Certainly Fibonacci ultra-thin magnetic films 
can be realized experimentally following the procedures 
of ref.[12] or [20] to grow the samples. 
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FIG. 1. The third and fifth Fibonacci generations and their 
magnetic counterpart. 
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FIG. 3. Magnetoresistance for the fifth (a) and seventh (b) 
Fibonacci generations with the same parameters of Fig. 2. In 
Fig. 3(a) the relative positions of magnetizations are indicated 
by the arrows, and the Fe double layer is indicated by the 
bigger arrow. 
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FIG. 2. Magnetization (a) and magnetoresistance (b) ver- 
sus magnetic field for the third Fibonacci generation with 
Hbi = —150 Oe and Htq = 50 Oe. In Fig. 2(a) the arrows 
indicate the relative positions of the magnetizations in each 
phase. 
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FIG. 4. Magnetoresistance for the third (a), fifth (b) and 
seventh (c) Fibonacci generations with Hm — —1.0 kOe and 
Htg = 0.1 kOe. 
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FIG. 5. Magnetoresistance for the third (a), fifth (b) and 
seventh (c) Fibonacci generations with \Hbi \ = Hi,g = 35 Oe, 
which correspond to a sample with Cr thickness equal to 25 
A. Note a striking self-similar pattern. In Fig. 5(a) and Fig. 
5(b) the arrows indicate the relative positions of the magne- 
tizations in each phase. 
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